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Chapter  1:  Introduction 

Over  the  past  decade  a  new  technology  has  emerged  using  surface 
acoustic  waves  (SAW)  for  signal  processing  applications  in  the  upper  VHP 
and  UHF  range.  These  applications  include  miniature  bandpass  filters, 
resonators,  convolvers,  correlators  and  Fourier  transformers  [  1  ].  The 
chief  advantage  of  acoustic  waves  over  their  electromagnetic  counterparts 
arises  from  the  fact  that  acoustic  velocities  in  solids  are  typically 
~  10  ^  times  the  velocity  of  light;  the  acoustic  wavelength  at  100  MHz  is 
about  30itn  compared  to  the  electromagnetic  wavelength  of  3  meters.  A 
significant  miniaturization  of  signal  processing  components  is  thus  achieved. 
Surface  waves  have  proved  particularly  useful  in  this  context  because  their 
energy  is  confined  within  a  wavelength  of  the  surface  and  is  accessible  all 
along  the  propagation  path.  Moreover,  unlike  most  guided  waves,  surface 
waves  are  non-dispersive . 

Most  signal  processing  devices  employ  surface  waves  in  piezo¬ 
electric  anisotropic  solids.  Periodic  arrays  of  electrodes  are  fabricated 
on  the  surface  to  interact  with  the  wave  and  produce  the  desired  results. 

An  accurate  description  of  these  interactions  is  essential  to  proper  filter 
design.  The  main  objective  of  this  thesis  is  to  develop,  from  basic  field 
theory,  simple  models  describing  these  interactions  that  can  be  used  by 
the  device  designer. 

The  following  chapters  deal  with  a  variety  of  interesting  acoustic 
field  problems.  Most  of  these  results  have  already  bean  published  and  the 
purpose  of  this  thesis  is  to  provide  a  compact  outline  of  the  approaches 
used.  Detailed  discussions  of  the  techniques  and  references  to  the  literature 
are  available  in  the  publications .  A  short  tabulation  of  the  contents  of  each 
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chapter  and  the  relevant  publications  are  listed  below: 

Chapter  Mo.  Description  Re  fere nee 

2  A  method  for  obtaining  the  velocity  and  field  2 

distributions  of  surface  waves  along  the  free 

surface  of  a  single  or  laminated  anisotropic, 
piezoelectric  medium. 

3  A  theoretical  and  experimental  study  of  non-  3,4, 5, 6 

dispersive  waves  along  an  edge. 

4  A  method  for  calculating  the  approximate  dis-  7,3 

persion  curves  for  surface  waves  in  a  medium 

with  periodic  surface  perturbations. 

5  A  simple  theory  based  on  the  concept  of  power  o 

exchange  for  describing  the  evolution  of  a 
propagating  wave  in  the  presence  of  source 

dis  tribucions . 

6  Derivation  of  the  piezoelectric  scattering  10 

matrix  for  surface  waves  from  an  electrode 

in  a  periodic  array, 

7  Derivation  of  the  mechanical  scattering  matrix  11,12 

for  surface  waves  from  an  electrode  in  a 

periodic  array. 

8  A  scatter  matrix  formulation  for  multistrip  13 

couplers  used  to  couple  surface  waves  from 

one  track  to  another. 

9  Transducer  structures  for  optimum  coupling  14 

to  surface  and  line  waves. 

10  Application  of  the  principle  of  superposition  15,16,17,13 

for  simplified  transducer  design. 


Chapters  2,  3  and  4  describe  a  non-iterative  variational  technique 
for  solving  the  acoustic  field  equations  as  applied  to  three  different 
boundary  conditions.  The  same  basic  technique  is  applied  in  Chapter  2  to 
describe  waves  confined  near  a  surface  and  ir.  Chapter  3  to  describe  waves 
confined  near  an  edge.  Chapter  4  describes  surface  waves  in  a  medium  with 
periodic  discontinuities  at  the  surface.  This  technique  has  bean  used  to 


calculate  the  dispersion  curves  for  surface  waves  in  a  periodic  array  of 
electrodes.  However,  it  is  sometimes  difficult  to  get  a  converged  solution 
using  a  reasonable  number  of  terms  in  the  expansion  of  the  field  quantities 
An  alternative  perturbation  approach  to  the  problem  is  described  in  the 
remainder  of  the  thesis. 

Chapter  5  describes  a  simple  perturbation  theory  that  is  used  in 
the  later  chapters  to  characterize  the  interaction  of  surface  waves  with 
a  periodic  electrode  array.  This  theory  describes  the  generation  of  surfac 
acoustic  waves  by  sources  such  as  charge  and  stress  distributions.  The 
scattering  problem  is  then  treated  as  a  process  of  generation  by  sources 
induced  by  the  interaction  of  an  unperturbed  wave  with  the  discontinuity. 

The  effect  of  an  electrode  in  a  periodic  array  on  a  propagating 
acoustic  wave  is  described  conveniently  in  terms  of  a  scatter  matrix  that 
relates  the  outgoing  acoustic  waves  and  the  current  into  the  electrode  to 
the  incoming  acoustic  waves  and  the  electrode  voltage  (Fig.  6.2,  Eq.  6.7). 
Chapters  6  and  7  evaluate  the  piezoelectric  and  the  mechanical  parts  of 
the  first  two  columns  (Eq.  6.7)  of  the  scatter  matrix  respectively,  while 
Chapters  9  and  10  involve  the  last  column  of  the  scatter  matrix.  Chapter  8 
describes  a  scatter  matrix  formulation  for  multistrip  couplers  which  are 
widely  used  in  surface  wave  devices  as  beam  compressors  and  track  changers. 
The  scatter  matrix  elements  are  deduced  from  the  results  of  Chapter  6. 
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Chapter  2:  Surface  Acoustic  Waves 

Surface  acoustic  waves  are  non-dispersivo  waves  along  the 
surface  of  a  solid  and  decay  exponentially  in  the  direction  normal  to  the 
surface  (Fig.  2.1).  The  velocity  and  field  distributions  for  surface  waves 
are  usually  obtained  using  an  iterative  procedure  to  satisfy  the  boundary 
conditions  at  the  interface  [  19  1 •  This  procedure  is  rather  complicated 
especially  in  layered  structures  [20  ]  where  boundary  conditions  have  to 
be  satisfied  at  several  interfaces  . 

This  chapter  describes  a  non-iterative  technique  for  obtaining 
surface  wave  solutions  in  anisotropic  piezoelectric  solids  [2  ].  The 

boundary  conditions  are  incorporated  into  the  acoustic  field  equations 
using  position-dependent  material  constants.  This  makes  it  possible  to 
account  for  the  boundary  conditions  directly  without  iteration  even  in 
laminated  media  with  several  interfaces.  To  find  the  velocity  and  field 
distributions  typically  requires  the  diagonalization  of  a  30x30  matrix 
which  takes  only  a  few  seconds  on  a  computer. 

The  coordinate  axes  x?  ,  x0  are  chosen  so  that  the  wave 

propagate^  in  the  x^-direction,  decays  in  the  x^-direction  and  is  uniform 
in  the  x^ -direc tion .  The  acoustic  field  is  described  by  the  particle 
displacements  in  the  three  directions  u^,  u,;  and  u^ .  In  piezoelectric 
solids  there  is  an  accompanying  electrostatic  field  described  by  the  elec¬ 
trostatic  potential,  2.  The  stress,  T  and  the  electrical  displacement,  D 
are  related  to  the  field  quantities  u  and  2  by  the  constitutive  relations: 
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where  C,  e  and  8  are  Che  stiffness,  piezoelectric  and  permittivity 
tensors  of  the  medium.  The  indices  a,  b,  c,  d  take  on  the  values  1,  2, 

3  and  summation  over  repeated  subscripts  is  implied.  The  field  equations 
are : 
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where  p  is  the  mass-density  of  the  solid.  The  first  equation  is  obtained 
by  equating  the  force  on  a  differential  element  to  its  mass  times  its 
acceleration  and  assuming  an  exp(j-jut)  dependence  for  all  field  quantities. 
The  second  equation  is  Laplace's  equation  from  electrostatics. 

Using  the  constitutive  relations  (2.1)  to  substitute  for  T  and  D 
in  the  field  equations  (2.2)  we  obtain  differential  equations  relating  the 
field  quantities  u  and  p).  For  a  homogeneous  medium  the  material  constants 
C,  e  and  8  are  constants  so  that  the  differential  equations  have  constant 
coefficients.  In  a  laminated  medium,  the  usual  procedure  is  to  solve 
these  equations  separately  tor  each  of  the  layers  and  then  look  for  the 
correct  wavenumber  at  which  the  normal  stresses  and  electrical  displacemen 
are  continuous  at  the  interfaces.  For  a  semi-infinite  solid  this  involves 
making  the  normal  stresses  zero  at  the  surface. 

In  the  present  approach  the  material  constants  are  written  as 
tunctions  of  x^  so  as  to  describe  the  composite  medium.  For  a  semi-infi¬ 
nite  solid  we  write, 


F(x2)  =  r 


•  > 


(2.3a') 
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where  F  denotes  any  of  the  material  constants  p,  C,  e  or  6  of  the  sub¬ 
strate.  F  is  its  position-dependent  value,  0(x2)  being  the  unit  step 
function.  For  a  layered  medium  we  write, 

?(x2)  =  F  9(x2-h)  +  F'[0(x2)  -  9(x2-h)]  (2.3b) 

where  F  is  the  material  constant  of  the  substrate  while  f'  is  that  of 
the  layer  on  top,  of  thickness  h. 

Using  this  method,  we  have  a  single  differential  equation  with 
variable  coefficients  for  the  composite  medium.  The  field  equations  in¬ 
volve  derivatives  of  the  material  constants  giving  rise  to  delta  functions 
at  the  interfaces  that  take  care  of  the  boundary  conditions  automatically 
[  21  ]. 

To  solve  this  composite  differential  equation,  the  field 
quantities  u  and  0  are  expanded  in  a  series  of  orthogonal  functions.  We 
know  that  the  excitations  are  wavelike  along  x^  and  decay  essentially  ex¬ 
ponentially  along  x2>  In  view  of  this,  the  field  quantities  are  expanded 
in  a  series  of  Laguerre  functions  in  x2>  The  Laguerre  functions  form  a 
complete  orthonormal  set  over  (O,®)  and  are  defined  as, 

|m(q)>  =  exp(-  .  Pm(q) /ml  (2.4) 

ch 

where  is  the  m  Laguerre  polynomial.  The  presence  of  the  exponen¬ 
tial  term  is  particularly  convenient  for  the  decaying  field  distributions 
we  are  seeking  to  describe.  It  is  found  that  five  terms  in  the  expansion 
provide  an  adequate  description  of  the  field.  In  the  propagation  direction 
the  fields  are  assumed  wavelike  with  a  exp(-jkx^)  dependence,  k  being  the 


wavenumber. 
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Using  this  expansion  the  differential  equation  is  converted 
to  an  algebraic  matrix  equation.  The  matrix  elements  are  determined  from 
integrals  involving  the  position-dependent  material  constants.  The  un¬ 
knowns  are  the  expansion  coefficients  of  the  various  field  quantities. 

The  lowest  eigenvalue  of  the  matrix  yields  the  velocity  of  the  surface 
wave  and  the  corresponding  eigenvector  gives  the  field  distributions. 

This  technique  has  been  applied  to  obtain  surface  wave  solutions 
for  various  orientations  of  lithium  niobate  (LiXbO., )  ;  dispersion  curves 
have  been  computed  for  practical  layered  structures  such  as  silicon 
dioxide  on  lithium  tnntalite  and  geld  on  LiUbO^ .  Fig.  2.2  shows  the 
surface  wave  velocity  against  k’n  where  k  is  the  wavenumber  and  h  is  the 
thickness  of  the  gold  layer.  This  is  in  agreement  with  results  reported 
earlier  using  iterative  techniques  [22].  The  details  have  been  described 
in  refs .  [2]  and  [ 7 ] . 


10 


Chapter  3:  Line  Acoustic  Waves 

The  guiding  of  surface  acoustic  waves  along  a  free  surface  involves 
only  one  boundary  in  contrast  to  most  wavaguiding  structures  which  require 
two  boundaries.  This  makes  surface  waves  non-dis persive  since  the  guiding 
structure  has  no  characteristic  dimension. 

The  guiding  action  of  the  surface  is  explained  qualitatively  in 
this  way.  The  particles  near  the  surface  have  greater  freedom  of  movement 
than  those  within  the  bulk  of  the  solid.  This  makes  the  surface  appear 
less  "stiff"  than  the  bulk  so  that  the  wavefront  moves  slower  near  the 
surface.  Since  waves  in  general  tend  to  be  guided  along  regions  of  lower 
velocity  the  surface  acts  as  a  guide. 

This  view  of  surface  guiding  leads  us  to  expect  that  there  will 
be  a  further  lowering  of  wave  velocity  near  an  edge  formed  by  two  stress- 
tree  surfaces .  These  waves,  called  line  acoustic  waves,  are  non-dispersive 
too  since  an  edge  has  no  characteristic  dimensions.  In  addition  line 
acoustic  waves  are  di f f ractiotiless  and  the  total  energy  remains  confined 
within  a  square  wavelength  from  the  edge. 


The  theoretical  analysis  of  surface  waves  described  in  chapter 
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lowest  velocity;  symmetric  nodes  where  the  two  tacos  move  in  and  out  in 
svnchronism  have  higher  velocities  and  in  most  cases  are  not  guided  modes 

Line  waves  are  characterized  by  high  power  density,  zero  dif¬ 
fraction  and  nondispersive  propagation  which  make  them  suitable  for  signa 
processing  applications .  However,  for  practical  device  applications  we 
need  edges  that  do  not  have  imperfections  any  larger  than  a  fraction  of  a 
wavelength;  in  the  URF  range  this  means  that  chips  and  striations  at  the 
edge  should  be  smaller  than  a  micron. 

Ws  were  able  to  obtain  edges  suitable  for  lossless  non-dispersi 
propagation  of  line  waves  up  to  at  least  220  MHz  (>.  ~  16 -am)  by  cleaving 
commercial  lithium  niobate  crystals  with  top  polished  surfaces  [4  ].  The 
cleavage  plane  in  LikbO^  (012)  lies  along  the  x-axis  at  an  angle  of  37.24 
to  the  y- axis  [2.3],  Using  different  cuts  of  LiNbO^  as  the  starting 
material,  wedges  with  different  interior  angles  are  obtained  as  shown  in 


Fig.  3.2. 

For  our  experiments  we  used  a  wedge  with  an  interior  angle  of 
84.9°  obtained  by  cleaving  127. SG'1  rotated  Y-cut  LiNbO^  (Fig.  3.3)  .  dine 
the  cleavage  plane  is  crystallojraphi cally  determined,  1 1 1 w  cl :i ^  1  c  or  r n g 
wedge  is  known  very  accurately.  The  theoretically  calculated  particle 
displacements  and  electrical  potentials  along  the  top  surface  are  shown 
in  Fig.  3 . --  for  a  100  mw  line  wave  along  this  wedge. 

Line  waves  at  40  MHz  were  excited  using  an  interdigital  trans¬ 
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The  measured  line  wave  velocity  was  3500  m/sec.  (+  17.),  which  compares  well 
with  the  theoretical  value  of  3488  m/sec. 

A  220  MHz  line  wave  device  [ 5  ]  was  built  using  single  phase 
transducers  (Fig.  3.6).  These  transducers  have  electrodes  A/2  wide  as  com¬ 
pared  to  the  X/8  wide  electrodes  required  for  reflection-free  interdigital 
transducers.  This  reduces  the  resolution  requirements  by  a  factor  of  4  so 
that  higher  frequency  devices  are  fabricated  with  the  same  photolithographic 

capability.  Fig.  3.7  shows  the  frequency  response  of  the  device  which  has 
2 

the  ideal  (Sinx/x)  shape  anticipated  for  a  non-dispersive  delay  line  with 
unweighted  transducers.  No  significant  losses  could  be  detected  along  a 
1  cm.  propagation  path. 

One  of  the  drawbacks  of  the  single  phase  transducer  is  its 
poorer  coupling  compared  to  the  interdigital  transducer.  A  sawtooth  trans¬ 
ducer  (Fig.  3.8)  combining  some  of  the  features  of  the  interdigital  and  the 
single  phase  transducer,  provides  better  coupling  [6  ].  A  discussion  of 
the  coupling  efficiency  of  various  types  of  transducers  is  contained  in 
chapter  9. 


Measured  frequency 
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Chapter  4:  SAW  in  Periodic  Laminated  Media 

In  chapter  2  an  analysis  of  surface  waves  in  laminated  media  was 
described  using  an  orthonormal  series  of  Laguerre  functions  to  express  field 
distributions  in  the  direction  perpendicular  to  the  surface.  However,  the 
medium  was  assumed  uniform  in  the  propagation  direction  x^,  so  that  the  field 
quantities  were  wavelike  in  x^  with  a  single  Fourier  component  exp(-jkx^) . 

We  have  extended  this  technique  to  include  periodic  variations  of  the  material 
constants  in  the  propagation  direction  (Fig.  4.1)  by  using  Fourier  series  to 
express  field  distributions  as  required  by  Floquet's  theorem: 

f(xl’*2>  =mVmn  ®m(x2)  expt~j(k+n  .  ^)  x^  (4.1) 

where  0  (x„)  is  the  Laguerre  function 
m  z 

p  =  length  of  one  period 

and  fmn  are  the  expansion  coefficients  for  the  field  quantity  f. 

As  before,  using  position-dependent  material  constants  to  describe  the 
composite  media  we  can  take  care  of  the  complex  boundary  conditions  in  this 
problem.  The  differential  equations  are  then  converted  to  matrix  equations 
and  the  phase  velocity  for  a  given  wavenumber  is  obtained  from  the  lowest 
matrix  eigenvalue  [  7 ] . 

The  chief  difficulty  of  this  method  lies  in  the  fact  that  many 
practical  problems  require  a  large  number  of  Fourier  terms  for  an  accurate 
description  of  the  fields.  This  makes  the  size  of  the  final  matrix  un¬ 
wieldy  for  the  numerical  computation  of  eigenvalues.  Perhaps  a  different 
basis  set  for  expressing  field  distributions  in  the  propagation  direction 
would  be  more  appropriate.  However,  this  has  not  been  investigated. 


The  technique  was  applied  to  various  practical  structures  with 
considerable  success.  Fig.  4.2  shows  the  velocity  for  different  wave- 
numbers,  k,  obtained  for  a  grating  of  thin  electrodes  on  Y-Z  LiNbO^ •  The 
wavenumber  is  normalized  to  j3  (=  2tt/p)  where  p  is  the  grating  period.  At 
k  =  .5Pq  there  is  a  stopband  as  anticipated.  The  velocity  is  discontinous 
at  this  point.  The  potential  and  charge  distributions  along  the  surface 
are  shown  in  Fig.  4.3  corresponding  to  the  two  values  of  velocity.  The 
flat  potential  regions  correspond  to  the  regions  covered  by  electrodes. 
Between  the  electrodes  the  charge  is  expected  to  be  zero.  As  can  be  seen 
from  the  figures,  this  boundary  condition  is  only  approximately  satisfied. 
This  is  because  only  a  finite  number  of  Fourier  terms  is  used.  The  present 
analysis  is  thus  a  variational  analysis  that  minimizes  the  velocity  eigen¬ 
value  for  an  assumed  expansion  of  field  quantities  in  orthogonal  functions. 

The  method  has  also  been  applied  to  grooved  arrays  and  thick 
electrodes.  These  results  are  described  in  [8  ].  This  technique  provides 
a  generally  applicable  scheme  for  calculating  approximate  dispersion  rela¬ 
tions  for  laminated  piezoelectric  structures  with  periodic  boundaries. 
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Figure  4.3:  Thin  conductive  grating  on  Y-Z  LiNbO^ 

(a)  Surface  potential  for  higher  velocity  solution 
at  k  =  .5  0o 

(b)  Surface  potential  for  lower  velocity  solution 
at  k  =  .5  0o 

(c)  Surface  charge  for  higher  velocity  solution 
at  k  =  .5  3o 

(d)  Surface  charge  for  lower  velocity  solution 
at  k  =  .5  0o 
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Chapter  5:  A  Simple  Perturbation  Theory  for  Propagating  Waves 

A  commonly  occurring  problem  in  wave  theory  is  to  describe  the 
changing  characteristics  of  a  propagating  wave  due  to  a  distribution  of 
sources  in  the  medium.  This  source  distribution  may  be  due  to  external 
generators,  as  in  transducers;  or  it  may  be  induced  by  the  incident  wave 
due  to  a  change  in  the  material  properties,  as  in  refleccors.  In  most 
practical  cases  the  perturbation  of  the  wave  is  small  over  a  wavelength, 
and  no  essential  change  takes  place  in  the  relationship  between  the  dif¬ 
ferent  field  components  comprising  the  mode.  Under  these  conditions  ever, 
a  complicated  coupled  mode  like  a  Rayleigh  wave  is  described  by  a  single 
complex  amplitude,  A.  An  unperturbed  wave  travelling  in  the  z -direction 
has  a  constant  amplitude  A;  in  the  presence  of  sources  the  wave  is  per¬ 
turbed  and  A  is  a  slowly  varying  function  of  z  described  by  a  simple 
differential  equation  (Fig.  5.1): 


]  dA  =  _1_  dP 
A  dz  2P  dz 


(5-1) 


a 

where  (1)  dp  is  the  power  (active  and/or  reactive)  flowing  from  the  source 
into  the  wave  in  a  distance  dz  along  the  propagation  path, 
and  (2)  ?  is  the  time -averaged  power  flow  of  the  propagating  wave. 

cl 

This  equation  follows  quite  easily  for  real  dP  from  the  relation 
2 

P  ~  ,a[  and  the  conservation  of  powder.  For  imaginary  (or  reactive)  dp 
the  relation  is  not  so  obvious  but  can  be  proved  by  considering  the  instan¬ 
taneous  (rather  than  time-averaged)  power  flow. 

Equation  (5.1)  is  simple  in  appearance  but  is  remarkably  versatil 
in  solving  the  type  of  perturbation  problems  we  have  been  discussing.  Its 
applicant . tty  ranges  from  the  simple  case  of  a  plane  wave  in  a  transmission 
line  to  the  “-ore  comp licated  case  of  guided  acoustic  waves  in  piezoelectric 
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media,  as  we  propose  Co  show  with  the  following  examples.  The  resu1 l3 
presented  here  are  not  new,  but  they  are  all  derived  from  the  same  equation 

(5.1) ,  thus  providing  a  simple  unified  physical  picture  for  a  wide  variety 
of  wave  perturbation  phenomena.  This  makes  it  straightforward  to  ex¬ 
tend  the  results  to  new  problems  involving  waves  with  more  complicated 

fie  ids  . 

1.  A  classical  transmission  line: 

We  start  with  a  trivial  example  for  which  exact  solutions  are 
available  from  well-known  conventional  techniques  and  show  how  Eq. 

(5.1)  predicts  the  same  results.  W'e  consider  a  transmission  line  modeled 
as  a  distributed  circuit  with  a  series  impedance,  Z  and  a  shunt  admit¬ 
tance,  Y.  The  characteristic  impedance,  Zq  and  the  propagation  constant. 


v  are  related  to  Z  and  Y  by, 
'0  J 


Zo  “  v  Y 


v  =  \ZY 
'o 


(5 . 2a) 


(5.2b) 


Now  if  we  change  the  shunt  admittance  by  an  amount  AY,  the  propagation 

constant  will  change  from  v  to  v  —  Av. 

J  o  'o 


Yq+  Ay  =  vZ  (Y  -  AY) 


AY, 


-  Yo  (1  +  2Y}  if  iY  <<  Y 


so  that, 


,  ’o  AY 

Ay  =  T  ‘  T 

provided  AY  is  much  less  than  Y;  that  is,  if  the  perturbation  is  sr 
Now  we  derive  Av  using  Zq.  (5.1).  The  additional  admit: 
AY  causes  a  shunt  current,  I  per  unit  length  given  by 


(5.3) 
ia  1 1 . 
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I  =  V  .  AY 
s 

where  V  is  the  potential  of  an  unperturbed  wave.  This  shunt  current  I 
may  be  considered  an  induced  "source"  and  the  power  flow  into  the  wave 
from  this  source  in  a  distance  dz  is  written  as, 


dP 


I  dz 
s 


AY  dz 


Using  dP  from  (5.4)  in  (5.1), 


_1  dA 
A  dz 


.  AY 


Integrating  (5.5),  we  have. 


A(z)  =  Aq  exp(-Ayz) 


where,  A  is  a  constant 
o 

and,  .  Zo  ._ 

Ay  =  -y  .  AY 


(5.4) 


(5.5) 


(5.6) 


(5.7) 


The  results  in  Eqs.  (5.3)  and  (5.7)  are  seen  to  be  equivalent,  noting 

from  Eq.  (5.2)  that  ZqY  =  Y  •  The  amplitude  A  used  here  can  be  any 

of  the  field  quantities  (voltage  or  current)  with  the  unperturbed  phase 

factor  exp(-y  z)  removed, 
o 

2.  Plane  bulk  acoustic  waves  in  a  conducting  medium: 

We  now  go  on  to  a  more  complicated  example.  Let  us  consider  a 
plane  bulk  acoustic  wave  propagating  in  an  insulating  solid  with  a  piezo- 
electrically  stiffened  velocity  v  .  The  wave  is  accompanied  by  a  traveling 

cL 

electrical  potential  $  .  We  wish  to  find  the  attenuation  and  the  change  in 
velocity  of  this  wave  if  the  solid  has  a  finite  conductivity  a  and  a  dc  field 
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is  impressed  t^.  Jrift  the  carriers  in  the  propagation  direction.  Exact 

solutions  for  this  situation  are  available  [25]  and  once  again  we  can 

compare  our  results  from  Eq.  (5.1)  with  well-known  results. 

In  a  conducting  medium  the  piezoelectric  potential  *  due  to 

a 

the  acoustic  wave  causes  the  charges  to  bunch,  forming  a  traveling  wave 

charge  distribution  p.  This  charge  distribution  produces  an  electrostatic 

potential,  3  that  is  obtained  from  Poisson's  eauation. 
r  e 

?  =  32  S  t  (5.8) 

o  e 

where  5  is  the  oernittivitv  of  the  medium  and  5  is  the  unperturbed  wave  - 

o 

number  of  the  acoustic  wave.  Tor  small  perturbations  we  may  ignore  the 

difference  between  3  and  the  perturbed  wavenumber  3  and  replace  '-/ ~z  bv 
'  o  ' 

-j  3  .  Since  we  are  considering  plane  waves,  we  have  set  Vix  =  p/ ~y  =  0. 

Let  J  be  the  current  per  unit  area  due  to  the  traveling  charge  d 
tribution.  This  current  is  considered  the  source  from  which  the  acoustic 
wave  absorbs  a  power,  dp,  per  unit  area  given  by, 


dt  2  -a  'bz-  • 


(5.9) 


From  current  continuity. 


or. 


CJ 


J  = 


-j  aij 


( 5 . 10a) 
(5. 10b) 


v  =  ■»!  3  being  the  acoustic  wave  velocity.  Using  (5.10a) 

a  '  o 


in  (5.9),  we  have. 


d-P 

dz 


4a  (jx?} 


(5.11) 


so  that ,  from  E  t .  '5.1), 
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1  dA  0a  ^ 

A  dz  4P„ 


(5.12) 


where  P  is  the  power  flow  per  unit  area  of  the  acoustic  wave.  In  this 

d. 

example  the  wave  in  the  insulating  solid  was  assumed  the  unperturbed  mode. 
The  acoustic  wave  thus  has  no  associated  charge  p  .  However,  in  general 
we  should  add  a  tern  (ju)Pa) *  0 to  the  numerator  in  (5.12). 

We  now  need  the  charge  distribution  p  in  terms  of  0a  to  complete 
the  analysis.  This  is  readily  done  by  writing  J  from  the  material  constants 
and  using  Eq.  (5.10b). 


J  =  v  p  =  j3  cr(0  +  0  )  +  v  p  +  j3  Dp 


(5.13) 


where, 


(1)  a  is  the  conductivity  of  the  medium 


(2)  D  is  the  diffusion  constant  of  the  carriers 

(3)  vd  is  the  drift  velocity  of  the  carriers  in  the  propagation 
direction  due  to  an  external  dc  field. 

Here  we  have  written  only  the  linear  component  of  current  having 
exp(-j3Qz)  dependence.  The  non-linear  acoustoelectric  component  involving 
the  product  of  0  and  p  is  neglected.  Using  Eq.  (5.8),  we  obtain 


from  (5.13), 


jrnp  = 


a0a 


-0  +  t) 


(5.14) 


where, 


y  =  1-  v,/v  (drift  parameter) 


(5.15a) 


id  =  ct/S  (conductivity  relaxation  frequency) 
2 

<d^=  va/D  (diffusion  frequency) 


(5.15b) 

(5.15c) 


Using  (5.14)  in  (5.12), 


*  ?) 
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so  that,  A  =  Aq  exp[-(a  +  j43)z] 


where 


a  .  -3 
+  J  a~ 
Ko  o 


<t  i 
■  a  a 


P  r- 

"oc 


4P 


Y-J 


(l.  15) 


X 

X 


Here  cz  is  the  attenuation  constant  and  Ap  is  the  change  in  the  wavenumber. 
Now,  for  a  perfectly  conducting  medium  with  no  drift  field,  we  can  obtain 


the  change  in  wavenumber  from  (5.16): 


43 


sc 


$  $ 
a  a 


4P 


3  x€ 
'  o 


(5.17) 


Using  (5.17)  in  (5.16), 


a 


,,  43 

.43  sc 

p  T  J  3  ~  3 

^o  ^o  'o 


x  /yw 
c 


J* 


1  - 


yjj 


(5.13) 


"cxd 


This  is  seen  to  be  the  same  result  as  in  [25’  noting  that  AS  /3  is  equal 


sc  '  o 


to  half  the  electromechanical  coupling  constant. 

3.  Surface  acoustic  waves  in  the  presence  of  distributed  sources: 

Now  we  come  to  an  example  for  which  no  exact  solutions  are 
available,  namely  the  propagation  of  a  surface  wave  in  the  presence  of 
stress  and  charge  distributions  on  the  surface  and  distributed  volume 
forces  and  charges  within  the  bulk. 

Consider  an  unperturbed  surface  wave  propagating  with  a  velocity 


cs  along  the  surface  of  a  semi-infinite  substrate  with  a  perfect  open 


circuit  ( 3  =  0)  above.  The  propagation  is  in  the  z-direction  and  the  wave 
is  assumed  uniform  in  the  x-direction.  Now  if  we  assume  stress  distribu¬ 


tions  T.  fz')  (i  =  x,v,z)  and  an  electrical  displacement  b  (z)  on  the  surfac 
ry  •  y 


(y  =  0)  then  the  power  d?^  flowing  into  the  surface  wave  per  unit  hear. 


33 


h 

i 

e 

X 


i 


width  in  a  distance  dz  is  given  by, 

dP  =  A*.  dz  [-  £  v*(y=0)  t'  +  0*(y=O)  jo®'  (5.19) 

b  ^  i=x,y,z  y 

where  A  is  the  amplitude  of  the  surface  wave  and  v  (y=0) ,  0(y=O)  repre¬ 
sent  the  particle  velocity  and  electrical  potential  at  the  surface  of  a 
unit  amplitude  unperturbed  surface  wave.  The  signs  of  the  terms  in  (5.19) 
correspond  to  the  y-direction  being  into  the  surface. 

If  there  is  a  body  force  density  F*  and  a  charge  density  p 
within  the  volume  then  there  is  an  additional  power  flow  dP^  into  the 
surface  wave  given  by, 

CD 

*  1 


dP 


V  -  1  •  2  j  dy  Li£c,y,z  \  K  +  «  (3«P')J 


(5.20) 


The  total  power  flow  into  the  wave  is. 


dP  =  dP  +  dP 
s  v 


(5.21) 


Now,  from  Eq.  (5.1)  we  have, 


K 


_L  ^A 

A  dz 


1  dP 

2  1 A 1 2  P  dZ 
1  n 


(5.22a) 


where  P^  is  the  power  carried  by  a  unit  amplitude  surface  wave,  so  that 
the  power  in  a  surface  wave  of  amplitude  A  is  given  by 


P  =  A  P  . 
a  1  1  n 


(5.22b) 


Using  (5.19),  (5.20)  and  (5.21)  in  (5.22a), 


4  p  —  =  f  +  f 
n  dz  s  v 


(5.23) 


where , 


f  =  -I  v .  (y=0)  T.'  +  fi’<(y=0)  jwD ^ 

s  i=x,y,z  iw  iy  w  J  y 


(5.24a) 
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f  =  dy  '  Z  v.  !•' '  -r  ;  Cjxz '  )  (5. 24b  i 

v  ..  _i=x,y,z  i  i  ' 

o  J 

Here  the  o rimed  quantises  represent  sources  while  unprined  quantities  are 

the  fields  due  to  the  unperturbed  wave.  As  noced  in  connection  with  bulk 

waves,  the  unperturbed  mode  is  assumed  not  to  have  any  stresses  or 

electrical  disolacement  D  at  the  surface;  nor  any  forces  F.  or  charges  p 

y  l 

within  the  volume.  Otherwise,  the  corresponding  power  flow  terms  (such 

as  ( j  xD  (v=0)  ;  ,  c  being  the  source  potential)  have  to  be  added  to 

e  anc  l  , 
s  v 

Sq.  (3.23)  is  really  the  well-known  normal -mode  equation  for 
surface  waves  that  was  derived  from  the  reciprocity  theorem  [26].  It  is 
the  starting  point  for  numerous  problems  involving  surface  wave  transducer 
reflectors  and  amplifiers.  It  is  interesting  to  note  that  it  follows  from 
the  trivial-looking  relation  (3.1).  The  appropriate  equation  for  more 
complicated  fields  as  in  guided  surface  waves  can  also  be  written  down  by 
inspection. 

To  summarize,  we  have  shown  that  the  perturbation  of  any  wave  by 
impressed  or  induced  sources  is  described  by  the  simple  relationship  state 
in  E'j.  (5.1).  The  application  of  this  principle  is  illustrated  -with  three 
examples-- the  classical  transmission  line,  the  propagation  of  plane  acoust 
waves  in  conductive  media  and  the  generation  of  SAW  by  impressed  stresses 
and  charges.  ■  believe  that  this  provides  a  simple  unified  picture  of 
wave  perturbation  that  will  prove  useful  in  the  treatment  of  newer  and 
more  complicated  problems . 
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Chapter  6:  Piezoelectric  Scatter  Matrix  for  Thin  Conductive  Gratings 

Periodic  arrays  of  thin  metallic  electrodes  are  widely  used  in 
surface  wave  devices  for  constructing  transducers  and  reflectors.  In 
reflectors  the  electrodes  are  either  left  unconnected  or  shorted  together 
to  form  a  periodic  grating  so  that  the  structure  is  truly  periodic.  In 
transducers,  however,  the  electrodes  are  connected  to  external  voltage 
sources  and  the  particular  sequence  and  phasing  of  electrode  voltages 
depends  on  the  desired  spectral  response  [27].  This  makes  the  structure 
electrically  non-periodic.  The  analysis  to  be  described  in  this  chapter 
assumes  electrical  periodicity  and  is  directly  applicable  to  reflectors 
and  multistrip  couplers.  A  modification  is  required  in  its  application 
to  transducers  which  is  presently  being  developed  [28], 

The  array  of  electrodes  at  the  surface  perturbs  the  boundary 
conditions  so  that  both  stress  and  charge  distributions  are  induced  at 
the  surface.  The  basic  problem  is  to  determine  these  induced  source  dis¬ 
tributions  due  to  an  incident  surface  wave;  the  perturbation  theory  developed 
in  the  last  chapter  (Eq.  5.23)  is  then  used  to  determine  the  phase  change 
of  the  incident  wave  and  the  generated  reflected  wave.  The  assumption, 
of  course,  is  that  the  perturbation  is  small  over  a  wavelength  which 
is  true  in  almost  all  useful  applications. 

In  this  chapter  we  will  consider  only  piezoelectric  scattering 
due  to  the  induced  charge  distribution;  the  mechanical  scattering  due  to 
the  induced  stress  distribution  is  discussed  separately  in  the  next  chapter. 

In  strong  piezoelectric  materials  like  LiNbO^  the  scattering  from  thin 
electrodes  is  predominantly  due  to  induced  charges;  but  in  weak  piezoelec¬ 
tric  materials  like  quartz  the  mechanical  scattering  is  more  important. 
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We  consider  .m  incident  surface  wave  from  the  left  on  a  periodic 
array  of  electrodes  each  of  which  is  connected  to  ground  through  an  imped¬ 
ance  Z  (Fig-  6.1).  The  electrical  potential  at  the  surface  associated 
with  the  wave  is  taken  as  a  measure  of  its  amplitude. 


s 


e-j  =  ,2 


(6.1) 


whe~e  3  is  the  unperturbed  wavenumber.  z  is  the  amplitude  ''comDlex)  of 
o  ■  +  ‘ 

the  wave  entering  an  electrode  from  the  left.  The  wave  leaving  to  the  right 
has  its  amplitude  changed  by  Am  .  In  addition  a  reflected  wave  of  amplitude 
As  is  generated  to  the  left.  At  and  As  are  obtained  by  integrating  Eq. 

-  -r 

5.23  over  one  period: 


Z  —  P  /  2 

O  "»  1  J  Z 

A*  -  ~ r-  D  (z)  eJ  •' 

+  2  -p/2  r 


where , 


As  = 


$  4> 

z  =  -t_L 

o  2P 


r-p/2 

V  -  i  Q  2 

v  Dv(z)  e 
-P/2  7 


(6.2a) 


(6.2b) 


?  being  the  power  carried  by  a  surface  acoustic  wave  with  a  surface 
potential  of  $  . 

To  determine  Dy(z) ,  we  assume  that  an  electrostatic  field  described 
by  Poisson's  equation  is  superposed  on  t'ne  piezoelectric  field.  The  potential 
<4  due  to  this  field  is  written  in  a  Fourier  series  from  Floquet's  theorem 
for  periodic  media: 

,(n) .  ,,(n). 


+  °> 


.(y»z) 


=  \ 


't  e 
L-i  n 
n=  -« 


! ^  ;  \y  P-jp 


(5.3) 


where. 


a(n) 


a  2tt 

3  4-  n. — 

o  p 


Reflection  and  transmission  of  surface  acoustic 
waves  at  an  electrode  in  a  periodic  array 
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The  exponential  decay  along  y  follows  from  Laplace' 
free  media. 

The  total  electrical  potential  at  the  sue 
-  i  3  z 

2(z)  =  3+e  J'o  +  -5e(y  =  0,  z) 


equation  for  charge- 


ace  is  given  by, 

(6.4) 


The  tangential  electric  field  at  the  surface  is  obtained  from  the  deriva¬ 
tive  of  the  potential: 


_  .  .  -j3  z  .  _(n)  -j3^z 

E7(z)  =  e  o  +  /  j3v  ’z  e 

n=  -  “ 


(6.5a) 


The  elec tros tatic  part  of  the  field  gives  rise  to  a  normal  electric  dis¬ 


placement  D  (z)  at  the  surface  that  is  written  from  Laplace's  eauation: 
y- 


D  (z)  =  (S  +  S  )  ) 

y  ^  p  o  L-, 


\  i "(n) i 


n= 


6 


-j3 


(n) 


(6.5b) 


where  is  the  effective  permittivity  of  the  substrate  and  is  the 

permittivity  of  the  medium  above  the  surface.  The  problem  now  is  to 

determine  the  a  's  such  that  (Fig.  6.1) 
n 


E  (z)  =  0 
z 


D  (z)  =  0 

y 


Id  <! 


I  <  M  <  p/2 


(6 . 6a) 
(6.6b) 


This  is  done  analytically  in  [10]  using  the  properties  of  the  Legendre 
polynomial.  Once  the  induced  charges  in  the  strips  D  (z)  is  determined 
bb  and  bb  are  readily  obtained  from  Ft.  (6.21  .  A  (3x3)  matrix 
relating  the  fields  at  the  two  acoustic  and  on-?  electrical  ports  (Fig. 
6.2)  is  obtained  analytically: 


f-  + 

A 

11 

A12 

A13 

k  - 
1  - 

As 

-u 

r  = 

A21 

A22 

A23 

l  I 

J 

__A31 

A32 

.\ 

“33_ 

(6.7) 
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Fig.  a. 2  is  labeled  in  terms  of  dimensionless  coordinates  s  and  9  in 

place  of  3  and  z.  In  the  above  discussion  we  assumed  a  wave,  t  incident 
o  + 

from  the  left.  The  second  column  of  the  A-matrix  corresponding  to  a  wave, 

3  incident  from  the  right  is  obtained  by  symmetry. 

The  (2x2)  scatter  matrix,  3,  relating  the  two  acouscic  ports  is 
easily  obtained  from  the  A-matrix  by  setting  I  -  v/z  where  z  is  the  ex¬ 
ternal  load  at  each  strip.  Unconnected  strips  are  described  by  Z  =  - 
while  shorted  strips  are  described  by  z  =  0. 

Thus  a  complete  scatter  matrix  describing  the  piezoelectric 
scattering  by  a  periodic  array  of  electrodes  is  obtained  analytically  at 
all  frequencies.  The  results  have  been  applied  to  determine  the  reflection 
coefficient  per  strip  at  different  stopbands  and  for  different  metaii- 
(=  a/p).  The  decails  are  described  in  [10]. 


zation  ratios 


41 


Chapter  7:  Mechanical  Scattering  from  Periodic  Arrays 

In  this  chapter  we  consider  the  scattering  due  to  a  thin  strip 

overlay  of  height  h  («  wavelength)  on  the  surface  (Fig.  7.1a).  An 

incident  Rayleigh  wave  along  a  free  surface  induces  normal  and  tangential 

stresses  TS  ,  TS  ,  TS  (Fig.  7.1b)  at  the  interface  between  the  strip  and 
xy  yy’  zy 

the  surface.  The  basic  problem  is  to  determine  these  stresses  in  terms 
of  the  particle  velocities  of  the  incident  wave;  the  stresses  can  then  be 
used  in  the  perturbation  Eq.  (5.23)  to  yield  the  phase  shift  of  the 
incident  wave  and  the  reflected  wave. 

To  a  first  approximation,  we  may  assume  that  the  strip  floats 
on  the  wave  so  that  the  entire  strip  moves  with  the  same  particle 
velocities  as  the  unperturbed  surface.  In  this  case  the  problem  amounts 
to  finding  the  forces  that  have  to  be  exerted  to  the  strip  by  the  inter¬ 
face  in  order  to  move  it  with  the  same  particle  velocities  v  ,  v  ,  v  as 

x  y  z 

the  unperturbed  surface.  This  is  done  by  considering  the  forces  acting 
on  a  differential  element  dxdz  of  the  strip  of  height  h.  From  the  x- 


directed  forces  (Fig.  7.2a), 


T  =  jtup  h  v 
xy  J  r  x 


(7.1a) 


surface  wave  particle  velocity  provided  we  can  express  T 


terms  of  v  ,  v  ,  v  .  This  is  done  aoproximata ly  from  the 
x  y  z 

equations  to  yield  [ 3 0 ] , 

'  1 
T  =  -  —  p.  v 
xz  C  x 

s 


1  4  U.  ( <\  +  p*  1 


plane  strain 


(7.2a) 

(7.2b) 


T'  =  0 
yz 


(7.2c) 


where  Cs  is  the  surface  wave  velocity  and  ,  a  are  the  Lame  constants 
of  the  strip  material.  For  anisotropic  scrips,  similar  equations  are 
derived  but  are  more  complicated. 


So  far  we  have  not  considered  the  effect  of  the  finite  strip 


width.  At  the 

edges  of  the  strip  (z  =  -  a/2) 

the 

situati 

on  changes 

so«r.e 

whac  from  that 

depicted  in  the  force  diagrams 

o  t 

Fig .  7.2 

.  There 

is 

no 

strip  material 

on  one  side  to  exert  the  stresses 

* 

2-  j 

m  ' 

X 

and  l' 

on 

the 

xz 

yz 

zz 

z-  face  of  the  strip.  So  in  order  to  keep  the  strip  vibrating  with  the 
same  particle  velocities,  the  substrate  has  to  exert  extra  forces  at 
the  edges  to  make  up  for  the  missing  material  to  the  left  at  z  =  -  a/2 
and  to  the  right  at  z  =  -  a/2. 


The  interface  stresses  thus  have  the  spatial  distribution 

shown  in  Fig.  7.3.  The  stresses  are  uniform  under  the  strip  with  delta 

,  '  .  4u,'0  \ 

functions  at  the  two  edges.  In  the  figure  a  =  u  and  a  =  - — -H - — — 

x  z  X  +  2u 

Using  these  stress  distributions  the  reflection  coefficients 
predicted  for  several  practical  cases  agree  well  with  measured  values.  The 
details  are  described  in  [11].  It  is  found,  however,  that  this  does  not 
predict  the  slight  lowering  of  phase  velocity  from  the  unperturbed  value 


that  is  experimentally  observed . 


ihe  physical  reason  for  this  velocity 


<P-1S8S 


Figure  7. 


{Ci 


Spatial  distribution  of 

(a)  x-directed  interface  stresses 

(b)  y-directed  interface  stresses 

(c)  z-directed  interface  stresses 


lowering  is  that ,  due  to  the  interface  stresses,  the  substrate  is  cl  is  tor 
somewhat,  and  the  assumption  of  the  strips  floating  on  the  wave  made  her 
is  not  strictly  true.  This  distortion  represents  a  relaxation  of  the 
system  causing  a  lowering  in  velocity.  Mathematically,  this  is  taken 
into  account  by  assuming  the  excitation  of  evanescent  modes  that  store 
energy' . 

The  first-order  stress  shown  in  Fig.  7.3  induces  particle  dis¬ 
placements  with  spatial  harmonics  separated  from  the  unperturbed  wave- 
number  p  by  an  integer  multiple  of  2~/p  where  p  is  the  period  (Floquet's 
Theorem)  of  the  array.  These  fields  store  energy  and  their  effect  on 
the  incident  wave  velocity  is  obtained  from  (5.23)  using  the  second-orde 
stress  generated  by  the  induced  particle  displacement  field.  A  complete 
analysis  of  this  scored  energy  effect  is  described  in  [12]  ;  the  phase 
velocity  and  reflection  are  correctly  predicted  using  this  analysis. 

The  analysis  described  in  this  chapter  yields  the  mechanical 
part  of  the  scatter  matrix  of  an  electrode  in  a  periodic  array.  These 
are  added  to  the  A^,  A  ,  A0 1  ,  A,  ?  obtained  in  chapter  5  (Eq.  6.7)  for 
the  piezoelectric  part  of  the  scattering.  The  other  matrix  elements  are 
however,  unchanged  by  mechanical  effects. 


47 


Chapter  8:  Multistrip  Couplers 

A  multistrip  coupler  (MSC)  is  a  periodic  array  of  electrodes 
that  is  used  to  couple  from  one  surface  wave  track  to  another  (Fig.  8.1). 
Each  strip  of  the  MSC  samples  the  potential  of  the  incoming  acoustic  wave 
and  applies  it  to  all  the  tracks  generating  acoustic  waves  in  each  track. 
Because  of  the  traveling  wave  nature  of  the  excitation,  the  generated  waves 
are  unidirectional;  reverse  waves  may  be  ignored  outside  the  stopbands. 
MSC's  are  widely  used  in  surface  wave  technology  as  track-changers,  beam 
compressors  and  couplers. 

In  this  chapter  we  will  derive  the  (MxM)  scatter  matrix  of  a 

single  strip  in  a  MSC  that  relates  the  incoming  and  outgoing  acoustic  waves 

in  the  M  tracks  (Fig.  8.2) 

M 


A?  =  )  b“  B“ 

i  k4i  ik  k 


,  i  =  1,  M 


(8.1) 


where  B^  is  the  incoming  wave  amplitude  in  the  kth  track  at  the  nth  strip 
and  A?  is  the  outgoing  wave  amplitude  in  the  ith  track  at  the  nth  strip. 
We  wish  to  determine  b^.  The  superscript  n  denoting  strip  number  is 
dropped  for  convenience.  The  scatter  matrix  of  successive  strips  can 
then  be  cascaded  to  obtain  the  scattering  characteristics  of  the  coupler 
as  a  whole. 


Neglecting  backward  waves,  we  may  write  in  terms  of  the  A- 
matrix  in  chapter  6, 


A  =  A(i)  B  +  A(i)  V 
i  A21  Bi  +  A23  V 


I  =  A*~^  B  +  A*^  V 
i  31  Hi  A33  V 


(8.2a) 


,  i  =  1,  M 


A 


(8,2b) 
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where  I.  is  Che  current  into  the  i  track,  V  is  the  common  voltage 

all  the  tracks,  and  A'  '  represents  the  A-matrix  elements  of  the  i*" 

mn 

If  a  load  Y  is  connected  externally  to  each  strip  then. 


+  VYl  =  0 


Using  I.  from  (8.2b)  in  (3.3), 


.1  A31}  Bi  +  V*  lYL+  . 

i=l  i=l 


So  that,  V  =  - 


A(i)  P 
A31  3i 


(8 .4a) 


where  Y  =  Y  +  }  A~^ 

Jj  iLu  3j 


(8.4b) 


Using  (3.4a)  in  (3.2a), 


a  =  A(i)  -  \  ( 1 )  \'  ,00  , 

i  21  Di  "  A23  ,  L  ,  A31  3k 

—  ^ 


(3.5) 


Comparing  w’ith  (3.1)  we  have  the  scatter  matrix  b  • 

X.  K. 

b  =  4(a>  5  -  A(i)  \(k) 

ik  A21  °ik  A23  A3I 


(S .  6) 


5  =  1 


i  =  k 


5  =  0 


The  first  term  in  (8.5)  denotes  the  direct  propagation  while  the  second 

term  denotes  coupling  from  other  tracks.  An  incident  wave  in  track  k 

(k5 

produces  a  current  proportional  to  A^' .  The  strip  voltage  is  related  to 
the  current  by  the  total  s trip  admittance  Y  .  The  strip  voltage  is  then  proportional 
to  A^jVy  which  generates  a  wave  in  track  i  proportional  toA^  ''Y. 
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A  simple  circuit  model  for  the  scatter  matrix  is  described  in  [13]. 
The  scatter  matrix  method  is  applicable  to  the  analysis  of  non-uniform 
couplers  with  slowly  varying  characteristics;  such  couplers  are  not  amen¬ 
able  to  analysis  by  conventional  techniques.  Details  of  this  technique 
with  examples  are  described  in  Reference  [13]. 
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Chapter  9:  A.  Uni  tied  Theory  of  Transducers 


In  the  last  three  chapters  we  have  discussed  the  scatter  matrix 
for  a  single  strip  in  a  periodic  array.  As  mentioned  earlier,  this  formu¬ 
lation  is  directly  applicable  to  reflectors  and  multistrip  couplers;  but 
modifications  are  necessary  to  apply  it  to  transducers  which  are  usually 
not  electrically  periodic.  A  complete  scatter  matrix  description  for 
transducers  is  presently  being  developed  [28]. 

The  scatter  matrix  approach  to  transducers,  though  accurate, 
requires  numerical  implementation  thus  lending  little  physical  insight. 

A  less  accurate  approach  known  as  the  "weak  coupling  approximation"  has 
the  advantage  of  yielding  simple  results  that  are  easily  interpreted 
physically.  The  acoustic  waves  generated  by  the  transducer  induce  charges 
in  the  electrodes  which  in  turn  regenerate  acoustic  waves.  This  regene¬ 
ration  is  neglected  in  the  weak-coupling  approximation.  With  this  approx¬ 
imation,  the  charge  distribution  in  the  array  is  calculated  from  electro¬ 
statics  alone.  The  problem  then  divides  into  two  steps: 

(1)  The  charge  distribution  on  the  electrodes  is  determined  from 
the  electrostatic  equation. 

(2)  The  calculated  charge  distribution  is  used  in  Eq.  (5.23) 
to  determine  the  generated  acoustic  waves. 

This  chapter  describes  the  second  step  in  terms  of  an  aribtrary  charge  dis¬ 
tribution.  The  next  chapter  discusses  the  first  step  for  SAW  interdigital 
transducers . 

In  this  chapter  we  will  discuss  the  radiation  conductance,  capa¬ 
citance  and  0  of  the  generalized  transducer  depicted  in  Fig.  9.1.  These 
results  are  universal  and  can  be  applied  to  the  excitation  of  any  acoustic 


mode  IS AW  or  LAW) ;  certain  unifying  features  common  to  all  transducers  are 


Figure  9.1:  A  generalized  transducer  producing  a  potential 

T  T 

distribution  V0  and  a  charge  distribution  VP 

proportional  to  the  terminal  voltage  V 
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revealed  through  this  general  formulation.  An  external  generator  impresses 
an  alternating  voltage  V  at  a  frequency  <u  across  the  terminals  of  the  trans¬ 
ducer  which  is  an  interconnected  pattern  of  electrodes  that  produces  a 

T 

potential  distribution  on  the  surface  given  by  Vs  (y,z).  The  function 
T 

3  (y,z)  is  characteristic  of  the  transducer  and  has  to  be  determined  from 
the  electrostatic  equations  as  discussed  earlier.  Corresponding  to  this 
potential  dis  tribution  we  have  a  charge  distribution  per  unit  area  in  the 
transducer  VpA(y,z).  The  peak  electrostatic  energy,  W  ,  stored  in  the 
transducer  is  given  by. 


wT  =  ~  <?Ti?T> 

where  \0T|,pT)  =  j’j  dydz  3T*(y,z)  p~(y,z) 
over  the 
transducer 


(9.1a) 


(9.1b) 


Since  the  peak  energy  in  a  capacitor  is  given  by  •jCV  ,  the  capacitance  of 
T 

the  transducer  C  is  obtained  as 


' ®  it  ’ 


(9.2) 


Mow,  let  us  obtain  the  radiation  conductance  of  the  transducers 
due  to  the  generation  of  any  particular  acoustic  mode.  Suppose  this  mode 
has  an  electrical  potential  at  the  surface  given  by  A3  (y,z)  where  A  is 
the  amplitude  of  the  wave.  The  total  power  carried  by  the  mode  is 
written  as  |a|  .P  .  We  then  have  from  equation  (5.23), 

Li 


*  L*  T 

dA  _  A  V  dv  $  (y,z)  ja;  p  (y,z) 

A  4 . 1a  1  2Pt 


(9.3) 


If  we  assume  the  wave  grows  from  zero  amplitude  to  a  finite  amplitude 


Aq  over  the  length  of  the  transducer,  then  is  obtained  by  integrating 


(9.3)  : 
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L  T 

where  (0  |p  )  is  defined  analogous  to  (9.1b). 
power,  P  ,  is  then  written  as, 

cl 


p 

a 


•  P 


L 


2  2 
cu  V 

16p7 


,  <0LlpT> I 


(9.4) 

The  radiated  acoustic 


(9.5) 


This  is  the  power  coupled  to  one  mode  in  one  direction  only.  There  will, 
in  general,  be  some  power  coupled  to  the  wave  in  the  other  direction,  as 
well  as  to  other  modes.  The  same  expression  (9.5)  will  hold  with  the 


appropriate  0  and  P  . 

Li 

Since  the  power  dissipated  in  a  conductance  G  is  given  by 

L 

1  2 

—  GjV  we  may  write  from  (9.5), 

CL=-§|-  I  <<6LlpT>I2  (9.6) 


Here  GT  represents  the  power  radiated  into  one  mode  in  one  direction 

Li  ' 1 11  "  ' " 

only.  The  expression  GT  may  be  written  in  terms  of  the  fractional  changes 
in  velocity  of  the  acoustic  mode  caused  by  completely  shorting  the  surface 
with  a  conducting  layer.  Suppose  that  a  charge  -A.pI'(y,z)  per  unit  area 
is  induced  on  the  conducting  layer  by  the  acoustic  wave  of  amplitude  A. 
Then  we  can  show  f r >m  (5.23)  that, 
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where  3^  Is  the  unperturbed  wavenumber  of  the  acoustic  mode.  Now  for  a 

travelling  wave  0^(y,z)  depends  on  z  through  an  exponential  phase  factor 

exp(-j3Qz) .  The  charge  induced  on  a  uniform  conducting  sheet  also  has  the 

L*  L 

same  phase  factor  so  that  0  (y,z)p  (y,z)  is  independent  of  z.  We  can  then 


wri te , 


L* 


j'dy  5>~  (y,z)  p“(y,z)  =  j-  (0“  1  ?  )  (9.8) 

where  L  is  the  length  of  the  transducer  and  (3  |p^)  is  defined  as  in 
(9.1b)  . 

Using  (9.7)  and  (9.8)  the  radiation  conductance  is  written  as, 
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We  now  have  the  radiation  conductance  and  capacitance  of  the  generalized 
transducer.  It  will  be  noted  here  that  if  more  than  one  acoustic  mode  is 
coupled  to  the  transducer  then  G  is  computed  similarly  for  each  mode  ail 

Lt 

the  conductances  are  placed  la  parallel.  Assuming  only  one  mode  we  have 
for  Q, 
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and  N  is  number  of  wavelengths  in  the  transducer  =  — 


(9.10) 


The  other  terms  are  redefined  here  for  convenience: 

(1)  7 ^ I p b  =  jY  dydz  2"(y,z)  o(y,z) 

over  the 
transducer 

T  T 

(2)  't  (y ,  z)  and  3  (y,z)  are  the  potential  and  charge  distributions 
produced  by  the  transducer  with  unit  voltage  across  its  terminals. 

(3)  2L(y,z)  is  the  potential  distribution  of  a  unit  amplitude  acoustic 


mode  along  a  free  surface. 
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variables  over  y,z  can  be  replaced  by  an  integration  of  the  Fourier  trans¬ 
forms  over  k  ,  k  ,  so  that  using  ("9.11b), 
v  z 
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(9.13a) 


where  •(  )  now  denotes  integration  over  k  ,  k  and  the  'oars  denote  resnec- 

y  z 

tive  Fourier  transforms.  The  two  terms  in  the  denominator  are  the  norms 
— T  — L  — 

of  p  and  o  with  a  weight  function  C  while  the  numerator  is  the  inner 

— L  —  t  — 

product  of  3  and  $  with  a  weight  function  C.  Clearly  F  is  maximum  when 

— L  —  T 

the  inner  product  is  maximum,  that  is,  when  3  =  3  .  This  is,  of  course, 
what  one  would  expect  intuitively.  We  note  here  that  F  can  be  written  in 


terms  of  charges  as  well: 

!(?L  C’liaT>!2 


<PT  c'V)  <pL  c'Mp1) 

The  quantity  F  is  a  good  figure  of  merit  for  a  transducer  exciting 
a  certain  acoustic  mode.  As  an  illustration  let  us  evaluate  it  for  a  SAW 
interdigitai  transducer  with  alternating  polarities.  Since  the  SAW  is  uni¬ 
form  along  y,  we  may  neglect  the  integrations  over  y  (or  k  ) .  Also  we  note 
that  for  fields  uniform  along  y, 

C(k  =  0,  k  )  =  6  <k  ] 
y  z  p 1  z 


(9.13b) 


(9. 14) 


where  S  is  the  effective  oermittivity  of  the  substrata.  This  is  shown 
P 

readily  from  Poisson's  equation. 

Suppose  the  transducer  potential  distribution  is  periodic  with  a 
Fourier  series  representation  to  the  form: 


3T(z)  =  Z  tT  e'j^O2 
n=  -=°  n 


(9.15) 


where  3  is  equal  to  the  wavenumber  of  the  acoustic  mode  being  excited. 
0 
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0L(z)  =  0L(O)  e'je0Z 


(9.16) 


Using  (9.14),  (9.15)  and  9.16)  we  obtain  C  from  (9.2),  G  from  (9.9) 

1  & 

and  F  from  (9.13a): 
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(9.17c) 


Here  the  Fourier  transform  versions  of  (9.2)  and  (9.9)  for  C  and  G 

T  3, 

have  been  used  and  the  summation  from  to  +°°  has  been  replaced  by 

— T  — T 

twice  the  summation  from  0  to  +®  assuming  0^  =  0  Also  G&  has  been 

multiplied  by  2  to  account  for  the  bidirectionality.  The  numerical 

value  of  F  was  computed  using  the  potential  distribution  for  a  SAW  IDT 

with  alternating  polarities  and  with  equal  electrodes  and  gaps.. 

It  is  evident  that  while  the  capacitance  arises  from  all  the 

T 

Fourier  components  of  the  transducer  distribution  function  0  ,  the  conduc¬ 
tance  comes  only  from  the  Fourier  component  that  matches  the  acoustic 

T 

wave  being  excited.  The  Q  is  thus  minimized  when  0  contains  only  one 

T  L 

Fourier  component  matching  the  acoustic  wave,  that  is,  when  0=0. 

The  equations  derived  in  this  chapter  are  applied  readily  to 
line  waves  as  well.  However,  for  any  given  transducer  we  have  to  obtain 
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Che  transducer  distribution  function  3  by  solving  the  electrostatic 
Poisson's  equation  with  the  appropriate  surface  boundary  conditions.  The 
line  wave  distribution  o'  is  obtained  from  the  field  equations  as  dis¬ 
cussed  in.  Chapter  3.  Fig.  3.4  shows  for  the  particular  wedge  used 
in  our  experiments. 

T 

The  problem  of  obtaining  3  for  the  line  wave  transducers  is  a 

three  dimensional  electrostatic  field  problem  and  unfortunately  this  has 

not  been  solved  for  the  single  phase  or  the  triangular  transducer  discussed 

in  Chapter  3.  However,  remembering  that  the  optimum  transducer  is  one  with 
L  T 

3  =  3  ,  we  can  get  a  feel  for  the  relative  Q's  of  different  types  of  trans¬ 
ducers  . 

T 

The  transducer  distribution  3  should  match  the  line  wave  dis¬ 
tribution  3^  both  along  the  propagation  direction  and  in  the  transverse 
direction.  Interdigital  transducers  provide  a  fairly  good  match  for  the 
propagation  direction  as  shown  by  the  high  F  computed  for  SAW  for  an 
interdigital  transducer.  However,  it  is  uniform  in  the  transverse  direction 
and  is  not  a  good  match  for  the  LAW  distribution  (Fig.  3.4).  In  fact, 
the  F  computed  for  LAW  for  a  1\  wide  IDT  is  ~  .15  showing  rather  poor 
efficiency.  The  sawtooth  transducer  is  similar  to  the  IDT  but  the  funda¬ 
mental  Fourier  component  is  weighted  across  the  beam  and  can  be  matched 
to  the  LAW  distribution  to  provide  a  higher  F.  However,  no  quantitative 
calculation  has  been  carried  out. 

The  single  phase  transducer  is  clearly  the  weakest  because  it 
provides  a  poor  match  even  in  the  propagation  direction.  Since  there  are 
no  negative  electrodes  between  the  positive  electrodes  the  potential  is 
not  forced  to  alternate  along  the  propagation  path  so  that  the  fundamental 
Fourier  component  is  rather  small. 


T  L 

The  condition  for  optimum  coupling  0=0  suggests  that  the 
transducer  in  Fig.  9.2  should  be  fairly  efficient.  However,  it  is 
difficult  to  implement  because  of  the  cross-overs  needed  to  make  con¬ 
nections  to  the  electrodes. 
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Chapter  10:  The  Superposition  Principle  as  Applied  to  Transducer  Analysis 
and  Design 

In  this  chapter  we  consider  an  interdigital  transducer  consisting 
of  a  periodic  array  of  electrodes  with  a  specified  voltage  sequence  (Fig. 
10.1).  As  mentioned  in  the  last  chapter,  the  generation  of  surface  waves 
by  the  transducers  is  computed  in  two  steps.  First,  the  charge  distri¬ 
bution  on  the  electrodes  is  determined  assuming  the  substrate  to  be  a  pure 
dielectric.  The  calculated  charge  distribution  is  then  used  in  Eq. 

(5.23)  to  determine  the  waves  generated  in  either  direction.  Here  we  will 
describe  an  application  of  the  superposition  principle  that  simplifies  the 
first  step  considerably  for  periodic  SAW  interdigital  transducers.  This 
first  step  of  the  problem  requires  us  to  solve  the  electrostatic  equations 
so  that  the  charge  is  zero  in  the  gaps  and  the  electrical  potential  has 
the  specified  values  on  the  electrodes.  This  is  a  non- trivial  problem 
especially  for  arbitrary  electrode  voltages.  Moreover,  the  problem  has 
to  be  solved  anew  for  each  new  set  of  voltages. 

In  this  chapter  we  will  show  that  the  charge  distribution  in  a 
periodic  transducer  is  written  as  a  convolution  product  of  the  electrode 
voltages  with  a  basic  charge  distribution  that  depends  only  on  the  elec¬ 
trode  metalization  ratio  (=a/p) ,  so  that  it  is  not  necessary  to  solve 
the  field  equations  for  each  new  set  of  voltages. 

Consider  the  transducer  in  Fig.  10.2  (a)  with  IV  applied  to  the 
central  electrode  and  all  other  electrodes  grounded.  The  charge  distri¬ 
bution  in  this  array  a(x)  (Fig.  10.2  (b))  is  the  basic  charge  distribution 
mentioned  earlier.  To  show  this  let  us  consider  the  three  cases  shown  in 
Fig.  10.3.  In  Fig.  10.3  (a)  the  charge  distribution  is  clearly  V^c(x) .  In 
Fig.  10.3  (b)  the  charge  distribution  is  V^a(x-p).  Now,  in  Fig.  10.3  (c)  the 
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Charge  distribution  in  a  periodic 
transducer  as  a  superposition  of 
basic  charge  distribution  functions 
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the  charge  distribution  is  obtained  by  superposing  the  charge  distributions 
in  Figs.  10.3  (a)  and  10.3  (b)  giving  VqO(x)  +  V^a(x-p) .  This  is  true  be¬ 
cause  the  field  equations  are  linear  and  the  sum  of  two  solutions  is  also  a 
valid  solution.  The  only  question  is  whether  it  has  the  correct  boundary 
values.  It  is  easy  to  check  that  it  does. 

Similarly  for  an  arbitrary  sequence  of  voltages  as  in  Fig.  10.1 
we  may  write  the  charge  distribution  Q(x)  as, 

Q(x)  =  2Vn  ct(x  -  np)  (10.1a) 

=  V(x)  *  o(x)  (10.1b) 

where  V(x)  =  £  6(x  -  np)  (10.1c) 

and  the  *  denotes  convolution.  It  will  be  noted  that  this  argument  is 
exact  only  for  an  infinite  array  of  electrodes  with  all  electrodes  occur¬ 
ring  periodically.  The  end  effects  are  easily  accounted  for  by  adding  a 
few  grounded  electrodes  at  each  end  since  the  effects  of  neighboring  elec¬ 
trodes  are  negligible  beyond  the  third  nearest  neighbor.  However,  if  any 
electrodes  are  withdrawn  from  the  array  (as  is  true  for  a  few  practical 
transducers)  then  this  analysis  is  not  applicable  in  its  present  form. 

It  can  be  shown  from  the  perturbation  equation  (5.23)  that  the 
total  surface  wave  amplitude,  Ac  at  a  frequency  uj  generated  by  a  distri- 
buted  charge  distribution  Q(x)  proportional  to  the  Fourier  transform  of 
the  charge  distribution  at  a  spatial  harmonic  equal  to  the  wavenumber  of 
the  surface  wave.  Thus, 

A  («)  «  <5(k  -  »/vg)  (10.2) 

where  Q(k)  is  the  Fourier  transform  of  Q(x) 
and  Vg  =  surface  wave  velocity. 
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From  (10.1)  we  than  have, 

Ag  (*)  «  V (k  =  aj/vs)  .  trfk  =  (10.3) 

where  V(k)  and  cr(k)  are  the  Fourier  transforms  of  V(x)  and  cr(x)  .  The 
overall  response  thus  splits  into  the  product  of  an  array  factor  V(k)  with 
an  element  factor  o(k) .  The  array  factor  is  obtained  from  a  simple  Fourier 
transformation  of  the  electrode  voltages  while  the  element  factor  is  ob¬ 
tained  once  and  for  all  from  field  theory.  Fig.  10.4  shows  the  element 
1/2  - 

factor  k  c(k)  plotted  against  k  for  various  values  of  the  metallization 
ratio.  These  plots  were  obtained  from  field  theory  [IS]. 

This  concept  has  been  applied  to  various  types  of  transducers 
and  the  details  are  described  in  [15],  [16],  and  [17].  Presently  the 
technique  is  being  extended  to  non-periodic  transducers  [29].  The 
principle  is  also  being  incorporated  into  a  general  scatter-matrix 
characterization  of  transducers  [28]. 


L 


70 


►  *- 

4 

H 

V 

i 


Re  ferences 


1.  L.  A.  Coldren  and  R.  L.  Rosenberg,  "SAW  Resonator  Filter  Overview: 

Design  and  Performance  Tradeoffs,"  IEEE  Ultrasonics  Symposium  1978. 

2.  5.  Datta  and  B.  J.  Hunsinger,  "Analysis  of  Surface  Waves  using 
Orthogonal  Functions,"  J.  Appl.  Phys . ,  Feb.  197S. 

3.  S.  Datta  and  3.  J.  Hunsinger,  "Analysis  of  Line  Acoustical  Waves  in 
General  Piezoelectric  Crystals,"  Phys.  Rev.  B  _15,  Nov.  1977. 

4.  S.  Datta,  M.  J.  Hoskins  and  B.  J.  Hunsinger,  "Line  Acoustic  Waves  on 
Cleaved  Edges,"  Appl.  Phvs.  Lett.,  1  Jan.  1978. 

5.  M.  J.  Hoskins,  5.  Datta  and  B.  J.  Hunsinger,  "UHF  Single-phase  Line 
Acoustic  Wave  Transducers,"  Appl.  Phys.  Lett.,  15  July  1973. 

6.  M.  J.  Hoskins,  5.  Datta,  C.  >1.  Panasik  and  3.  J.  Hunsinger,  "Line 

Acoustic  Waves  on  Cleaved  Edges  in  LiNbO,  and  GaAs , "  IEEE  Ultrasonics 
Symposium  1973.  ' 

7.  S.  Datta  and  3.  J.  Hunsinger,  "Surface  Acoustic  Waves  in  Periodic 
Structures,"  J.  Appl.  Phys.,  to  be  published. 

8.  S.  Datta  and  B.  J.  Hunsinger,  "Reflection  of  Surface  Acoustic  Waves  from 
Periodic  Structures  in  Layered  Media,"  IEEE  Ultrasonics  Symposium  197S. 

9.  S.  Datta  and  3.  J.  Hunsinger,  "A  Simple  Perturbation  Theory  for 
Propagating  Waves,"  submitted  to  J.  Appl.  Phys. 

10.  S.  Datta  and  3.  J.  Hunsinger,  "An  Analytical  Theory  for  Piezoelectric 
Scattering  by  Periodic  Arrays,"  submitted  to  J.  Appl.  Phys. 

11.  S.  Datta  and  3.  J.  Hunsinger,  "First-order  Reflection  Coefficient  of 
Surface  Acoustic  Waves  from  Thin  Strip  Overlays,"  J.  Appl.  Phys., 

to  be  published. 

12.  S.  Datta  and  B.  J.  Hunsinger,  "An  Analysis  of  Energy  Storage  Effects 
on  SAW  Propagation  in  Periodic  Arrays,"  submitted  to  IEEE  Trans,  on 
Sonics  and  Ultrasonics. 

13.  S.  Datta  and  B.  J.  Hunsinger,  "A  Model  for  Multitrack  Non-periodic 
Multistrip  Couplers,"  J.  Appl.  Phys.,  July  1978. 

14.  S.  Datta  and  B.  J.  Hunsinger,  "A  Unified  Theory  of  Transducers," 
unpublished . 


15.  S.  Datta  and  3.  J.  Hunsinger,  "Redefined  Element  Factor  for  Simplified 
IDT  Design,"  Electron.  Lett.,  9  Nov.  1978. 


71 


16.  B.  J.  Hunsinger  and  S.  Datta,  "A  Generalized  Model  for  Periodic  Trans¬ 
ducers  with  Arbitrary  Voltages,"  IEEE  Ultrasonics  Symposium  1978. 

17.  S.  Datta,  B.  J.  Hunsinger,  and  D.  C.  Malocha,  "A  Generalized  Model  for 
Periodic  Transducers  with  Arbitrary  Voltages,"  IEEE  Trans,  on  Sonics 
and  Ultrasonics,  May  1979. 

18.  S.  Datta  and  B.  J.  Hunsinger,  "Element  Factor  for  Periodic  Transducers," 
IEEE  Trans,  on  Sonics  and  Ultrasonics)  to  be  published. 

19.  W.  R.  Jones,  J.  J.  Campbell  and  S.  L.  Veilleux,  Final  Report  to  AFCRL, 
Contract  No.  F19628-69-C-0132 . 

20.  L.  P.  Solie,  "Piezoelectric  Waves  on  Layered  Substrates,"  J.  Appl.  Phys . , 
February  1973. 

21.  A.  A.  Maradudin,  E.  W.  Montroll,  G.  H.  Weiss,  I.  P.  Ipatova,  Theory  of 
Lattice  Dynamics  in  the  Harmonic  Approximation,  2nd  ed.  (Academic  Press, 
New  York,  1971),  p.  524. 

22.  B.  A.  Auld,  Acoustic  Fields  and  Waves  in  Solids  (Wiley,  New  York,  1973), 
Vol.  II,  p.  282. 

23.  W.  L.  Bond,  Crystal  Technology  (Wiley,  New  York,  1976). 

24.  R.  L.  Whitman  and  A.  Korpel,  "Probing  of  Acoustic  Surface  Perturbations 
by  Coherent  Light,"  Appl.  Optics,  August  1969. 

25.  J.  H.  McFee,  "Transmission  and  Amplification  of  Acoustic  Waves  in  Piezo¬ 
electric  Semiconductors,"  Physical  Acoustics,  Vol.  IV,  Part  A,  ed.  by 

W.  P.  Mason  (Academic  Press,  New  York  and  London,  1966). 

26.  B.  A.  Auld,  Acoustic  Fields  and  Waves  in  Solids,  Vol.  II  (Wiley  &  Sons, 
Inc.,  1973),  p.  162. 

27.  C.  S.  Hartmann,  D.  T.  Bell  and  R,  S.  Rosenfeld,  "Impulse  Model  Design 
of  Acoustic  Surface-Wave  Filters,"  IEEE  Trans.  MTT-21  (1973). 

28.  C.  M.  Panasik  and  B.  J.  Hunsinger,  "Scatter  Matrix-based  Transducer 
Model,"  submitted  to  IEEE  Trans,  on  Sonics  and  Ultrasonics. 

29.  A.  L.  Lentine,  S.  Datta  and  B.  J.  Hunsinger,  "Charge  Distribution  for 
Non-periodic  Transducers  Using  a  Circuit  Model,"  submitted  to  IEEE 
Trans,  on  Sonics  and  Ultrasonics. 

30.  B.  A.  Auld,  Acoustic  Fields  and  Waves  in  Solids  (Wiley,  New  York,  1973), 
Vol.  II,  p.  277. 


72 


Vita 

Supriyo  Datta  was  born  in  Dibrugarh,  India  in  1954.  He  obtained 
his  5. Tech,  degree  from  the  Indian  Institute  of  Technology,  Kharagpur  in 
1975  and  his  M.S.  degree  from  the  University  of  Illinois  at  Urbana- Champaign 
in  1977. 


